Induced Magnetic Ordering by Proton Irradiation in Graphite 
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We provide evidence that proton irradiation of energy 2.25 MeV on highly-oriented pyrolytic 
graphite samples triggers ferro- or ferrimagnetism. Measurements performed with a superconducting 
quantum interferometer device (SQUID) and magnetic force microscopy (MFM) reveal that the 
magnetic ordering is stable at room temperature. 



The search for macroscopic magnetic ordering phenom- 
ena in organic materials has been pushed to the fron- 
tiers of physics, chemistry and materials science since 
the discovery of ferromagnetism in the open-shell rad- 
ical p-nitrophenylnitronyl nitroxide (p— NPNN) |1| and 
in tetrakis(dimethylamino)ethylene [TDAE] + Cgo charge- 
transfer salt [2j with Curie temperatures T c = 0.6 K 
and 16 K, respectively. Although a remarkable amount 
of studies have been performed in the last years build- 
ing molecules with unpaired 7r-electrons |3] , the scientific 
community did not become aware of the possible exis- 
tence of a room temperature metal- free organic ferromag- 
net till the accidental "discovery" of ferromagnetism in 
polymerized fullerenes was reported 0,0 ■ The study of 
the magnetic properties of polymerized fullcrcnc was trig- 
gered by the previously reported weak ferromagnetic-like 
magnetization loops in highly-oriented pyrolytic graphite 
(HOPG) [3 that indicate the existence of magnetic order- 
ing with T c much above room temperature [jj . The origin 
for the reported ferro- or ferrimagnetism in polymerized 
Ceo as well as in HOPG remains unclear. The contri- 
bution of ferromagnetic impurities, like iron, has been 
thoroughly investigated @>13 and it is something one 
should carefully take into account, specially when the fer- 
romagnetic signal under study is small. Ferromagnetism 
in a graphite nodule from a meteorite with T c ~ 570 K 
has been recently suggested to be due to a large mag- 
netic proximity effect at the interface with magnetite or 
kamacite inclusions As in Ref. @ most of the mag- 
netism specialists doubt that a room-temperature mag- 
net can exist in a material with only s- and p-electrons. 
Clear evidence for its existence and reproducibility as 
well as the report of a relatively easy method to pro- 
duce it, is of fundamental interest for a broad spectrum 
of natural sciences. 

Early reports on room-temperature ferromagnetic be- 
havior in some carbon-based structures were apparently 
not taken seriously by the scientific community, partially 
due to the weakness of the ferromagnetic signals added 
to the unknown contribution of magnetic impurities and 
low reproducibility [ldj . From those early works our at- 
tention was focused to the magnetic properties found in 
amorphous-like carbon prepared from different hydrogen- 
rich starting materials where an increase of the saturation 
magnetization with the hy drog en concentration in the 
starting material was found jllj . The origin for the mag- 



netic ordering has been related to the mixture of carbon 
atoms with sp 2 and sp 3 bonds, which was pre dicted to 
reach a magnetization higher than in pure Fe |l2j • Hydro- 
gen, on the other hand, was assumed to have a role only 
in the formation of the amorphous carbon structure . 
New theoretical predictions, however, show that hydro- 
genated graphite can display spontaneous magnetization 
coming from different numbers of mono- and dihydro- 
genated carbon atoms We show in this letter that 
implanted protons in HOPG triggers ferromagnetic (or 
fcrrimagnetic) ordering with a Curie temperature above 
room temperature. 

Two HOPG (ZYA grade) samples from Advanced Ce- 
ramics Co. with a content of magnetic metallic impurities 
below 1 ppm were irradiated with a 2.25 MeV proton 
microbeam. The samples' dimension was 2 x 2 x 0.1 mm 3 . 
Before irradiation they were glued with varnish on a 
high-purity Si substrate and the magnetic moment of 
the whole ensemble as well as of the Si substrate alone 
were measured. Due to the small sample mass the con- 
tribution of the HOPG virgin sample ferromagnetic sig- 
nal (saturation magnetization M s ~ 3 x 10~ 4 emu/g 
at temperatures T < 300 K 0) is overwhelmed by the 
diamagnetic signal of the Si-substrate for fields applied 
parallel to the graphene layers. Rutherford Backscat- 
tering Spectroscopy and Particle Induced X-ray Emis- 
sion spectra were recorded simultaneously with the ir- 
radiation allowing us to check the purity of the sam- 
ple at the different irradiation stages. Assuming that 
all the (measured) Fe impurities (< 0.3/ig Fe per gram 
graphite) in our samples behave as bulk ferromagnetic 
material, a rather unrealistic assumption, their maxi- 
mum contribution to the magnetic moment in our sam- 
ples would be to < 6.1 x 10 -8 emu. The irradiated and 
non-irradiated areas were further characterized by atomic 
force and magnetic force microscopy (AFM,MFM) with 
a Nanoscope III scanning probe microscope operating 
in "tapping/lift" mode, using standard tips coated with 
a magnetic CoCr film magnetized normal to the sam- 
ple surface. The magnetization measurements were per- 
formed with a SQUID magnetometer from Quantum De- 
sign with the reciprocating sample option (RSO) and a 
sensitivity of < 10 _7 emu. 

Four irradiations were consecutively applied to sample 
1, namely: stage #1: homogeneous irradiation of an area 
1720 x 1720 ^m 2 , dose: 0.99 pC/^m 2 , total charge: 2.93 
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fiC; #2: 100 x 100 spots of 2/zm diameter each, on an 
area 570 x 570 /zm 2 in the middle of the sample, dose: 0.3 
nC//Ltm 2 , total charge: ~ 8 /j,C; #3: 4 spots of 0.8 mm di- 
ameter each, dose: 0.3 nC//mi 2 , total charge: ~ 600 /uG; 
#4: the same as #3. To check the reproducibility of 
our procedure as well as to rule out possible contamina- 
tion during the handling of the sample, a new piece of 
the virgin HOPG sample was prepared in a similar way 
and fixed to a different Si-substrate. On this sample 2 
three spots of 0.8 mm diameter, dose: 0.3 nC//im 2 , total 
charge: 450 [iC were irradiated, similar to stage #3 but 
one spot less. 

To estimate the defect density created by the proton 
beam, Monte Carlo simulations (SRIM2003 14]) were 
performed with full damage cascades assuming a dis- 
placement energy of Ed = 35 eV for the creation of 
a Frenkel pair in graphite. According to these calcu- 
lations, 2.25 MeV protons come to rest at a depth of 
~ 46 jum from the graphite surface with a 1.75 fxm width 
(full width at half maximum of the implantation profile) . 
Since hydrogen is fully trapped in HOPG up to a con- 
centration of 0.45 H/C [l5j, all implanted protons will be 
trapped for doses up to 14 nC//im 2 when they come to 
rest in the HOPG sample. Assuming that no annealing of 
defects during irradiation occurs and that the damaged 
regions are of the size of the irradiated one, we obtain 
that for proton dose of 0.3 nC//zm 2 the defect density is 
~ 2 x 10 20 cm 3 , i.e. only ~ 0.1% of the carbon atoms are 
displaced. For a dose 250 times higher we expect that 
35% of the carbon atoms are displaced and the carbon 
lattice is amorphous. 

Figure ^a) shows the measured magnetic moment of 
sample 1 after irradiation #1, which is similar to the vir- 
gin sample within the scale of the figure. The field was 
applied parallel to the graphene planes. The main part 
(~ 90%) of the diamagnetic signal is due to the Si sub- 
strate. In the same figure we show the magnetic moment 
of the same ensemble after irradiation #3 where we can 
clearly recognize the s-shape curve, see Fig. ^a), with- 
out any background substraction. After subtraction of 
the magnetic moment of the substrate we obtain for sam- 
ple 1 the results depicted in Fig. ^b) . We observe a clear 
increase of the ferromagnetic loop with irradiation steps 
#2 and #3. For the last step #4 the magnetic moment 
of the irradiated sample decreased slightly. The coercive 
fields are between 80 and 110 Oe and weak tempera- 
ture dependent. The remanent magnetic moments are 
approximately 20% of the values of the saturation mag- 
netic moments. For the sample with the highest magnetic 
moment (sample 1 after stage #3): m rcm ~ 1 x 10~ 6 emu. 

The observed changes of the magnetic moment after 
the different irradiation stages are due to the irradiation 
and not to different misalignments of the sample position 
respect to the applied field. To check this we measured 
sample 1 in the irradiation stage #3 and #4 for the other 
field direction (parallel to the c-axis of graphite). After 
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FIG. 1: The magnetic moment (in units of 10 6 emu = 
10~ 9 Am 2 ) measured at T — 300K as a function of magnetic 
field (IkOe = 10 6 /47rAm _1 ) by cycling the field from zero to 
+10 kOe, from +10 kOe to -10 kOe and back to +10 kOe 
for sample 1 glued on a silicon substrate, after various proton 
irradiations, (a) Total magnetic moment without any back- 
ground subtraction for the sample after stage #1 (□) and 
#3 (•) irradiations, (b) Magnetic moment after subtraction 
of the sample holder magnetic moment, for the sample after 
the first (□), second (*), third (•) and fourth (v) irradiation 
stages. 



subtraction of the diamagnetic signal from graphite and 
Si substrate, the magnetic moment loops were similar as 
for the other field direction; this is an indication for a 
low anisotropy of the magnetic ordering produced by the 
irradiation. 

Figure Ufa) shows the magnetic moment of sample 
2 before and after irradiation of the three spots. The 
background signal due to the substrate was subtracted. 
We see a clear development of the ferromagnetic signal 
after irradiation. The saturation magnetic moment is 
slightly smaller than the one we would expect assuming 
a linear relation between saturation moment m s and to- 
tal implanted charge C t , namely ~ (3/4) x 5 x 10~ 6 = 
3.8 x 10 -6 emu. This value is 10% larger than the mea- 
sured value m s ~ 3.5 x 10~ 6 emu. Although the num- 
ber of points is too small to provide with certainty the 
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FIG. 2: (a) Magnetic moment (in units of 10 _6 emu) mea- 
sured at T = 300K as a function of magnetic field as in 
Fig. 1(b), for sample 2 before (□) and after (•) proton ir- 
radiation, (b) Measured saturation magnetic moment as 
a function of the total irradiated charge Ct for sample 1 
(•) and for sample 2 (□). The dashed line is the funtion 



= 0.22[10- 6 emu//iC - 5 ]C t 



relationship between m s and Ct, our results indicate a 
relation of the type m s oc y/Ut, see Fig.|2JT>). It is inter- 
esting to note that after the #4 irradiation the magnetic 
moment of sample 1 decreases, see Fig. Ufa), indicating 
that there might be a competition between the produced 
disorder and the implanted charge, which determines the 
total magnetic ordering. 

Hysteresis loops were measured at 5 K, 300 K and 
380 K. In this temperature range there is no signifi- 
cant change of the ferromagnetic loops with temperature. 
This result is similar to that obtained for the "intrinsic'' 
ferromagnetic-like loops of graphite in Ref. Q , and indi- 
cates a Curie temperature much above 400 K, which is 
of importance for future applications. 

The MFM pictures obtained at the irradiated areas 
show homogeneously distributed magnetic domains, as 
depicted in Fig. |21 The periodicity and/or the width 
of the domains depends on the irradiation dose, rang- 
ing from 4 fj,m in the second irradiated stage area to 
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FIG. 3: Top: Phase gradient pictures obtained from MFM 
at three surfaces of sample 1 corresponding to the irradiation 
stages #1, #2 and #3, see sketch at the bottom left of the 
figure. Bottom right: Topography and phase gradient line 
scans of the middle MFM picture obtained at stage #2. 



about 2/im at the area of the third stage. A weak do- 
main structure is recognized at the area of stage #1, see 
Fig. |21 We stress that the magnetic signal is not cor- 
related to the topography in any irradiation stage. For 
comparison, in Fig. |31 we show a topography and mag- 
netic line-scans obtained normal to the magnetic domain 
structure at stage #2. The MFM pictures shown in Fig-El 
were obtained at a distance of 50 nm between tip and 
surface. Due to the small coercive field of the magnetic 
surface and the influence of the magnetic tip we have 
observed that the magnetic domain distribution depends 
on the distance between tip and surface. Theoretically, 
the phase shift A<j> measured by MFM should be pro- 
portional to the magnetic force gradient, which should 
depend on the distance between the tip apex and sam- 
ple surface h as (h + 5)~ 5 , being 8 the distance between 
tip apex and the effective position of the magnetic mo- 
ment in the tip. Measuring Ac/> as a function of h we 
estimate S ~ (900 ± 100) nm. With this value we esti- 
mate the maximum magnetic moment observed by the 
tip m~3x 10~ 15 Am 2 at the area of largest magnetic 
moment (stages #3 and #4). We note that measure- 
ments on different magnetic samples indicate that S de- 
pends on the magnetic properties of the measured sur- 
face and therefore the calculated magnetic moment from 
MFM data should be taken as a rough estimate only 
[lfij . Assuming that the magnetic moment measured by 
MFM corresponds to a region of area ~ 10~ 12 m 2 and 
~ 1/nm depth, we can estimate the total saturation mag- 
netic moment m s the SQUID would measure at stage #3 
asm s - [3xl0" 15 x4x7rx(0.4xl0" 3 ) 2 /10" 12 ]10 3 emu= 
6 x 10 -6 emu, in reasonable agreement with the SQUID 
measurements. 

Although a proton penetration depth of ~ 46 /nm is 
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expected, we do not know how large is the depth for the 
ferromagnetic active part of the irradiated sample. This 
will be studied by irradiating samples of different thick- 
ness in the future. If we assume a depth of 1 ^tm we 
get a magnetization M ~ 1.1 emu/g for sample 1 after 
stage #3, which is ~ 1% of that of magnetite and of the 
same order as that obtained for amorphous-like carbon 
prepared by CVD from hydrogen-rich targets yjj . From 
MFM, the maximum phase shift observed in the irradi- 
ated regions is of the order of A<^ max ~ 0.3°. From mea- 
surements on regions with Fe impurities (mainly Fe3 04 
of micrometer size) with the same tip and arrangement 
[lfij ] we obtain A0 max ~ 10° . . . 20°. That means that 
A^max obtained in the irradiated regions #3 is of the 
order of 1% of that of magnetite, in agreement with the 
rough estimate made above using the SQUID data and 
assuming 1 /im ferromagnetic active thickness. 

Magnetic force microscopy measurements on proton- 
irradiated spots on HOPG show that it is possible to cre- 
ate ferromagnetic active areas of ~ 1 /im 2 with very sharp 
borders of the magnetic force gradient, created mainly at 
low proton doses [l7| . That means that it is possible to 
produce magnetic microstructures of arbitrary shape by 
a dedicated proton beam scanning system. To check the 
importance of hydrogen, we have produced spots of sim- 
ilar topography with 1.5 MeV helium ion irradiation at 
different doses (> 0.05 nC//im 2 ). Those irradiated re- 
gions showed no significant magnetic signal in MFM and 
stress the significance of hydro gen in the formation of the 
magnetic ordering in graphite |17| . 

From SQUID and MFM measurements we conclude 
that magnetic ordering appears after proton irradiation 
in graphite and that this is stable at room temperature. 
Our results provide a possible answer to the origin of fer- 
romagnetism in other carbon-based structures 0, 0, , 
which is not correlated with magnetic impurities. Future 
experiments should clarify the lattice position of hydro- 
gen in graphite as well as its role in the magnetic or- 
dering. The characterization of our samples with Micro- 
Raman as well as XPS is of particular interest and will 
be performed in the future. Systematic investigations of 
the magnetic properties of different carbon materials ir- 
radiated at different ion doses are necessary to gain more 
insight into the physical processes involved as well as to 
broader the spectrum for future applications. MFM mea- 
surements on proton irradiated amorphous-carbon films 
indicate the existence of weak magnetic ordering. As 
possible explanation for our results we refer to the work 
in Ref. |l3j where a sp 3 -sp 2 ferrimagnetic structure can 
arise in a graphene plane through the existence of mono- 
and dihydrogenated carbon atoms. Future work should 



provide an answer to the question to what extent the 
graphite structure represents a quantitative advantage 
for the formation of the magnetic ordering. 

Our experimental findings open up a new research field 
in magnetism with possible applications in spin electron- 
ics, since according to Ref. |l3( the electrons in the hy- 
drogenated graphene should have a very large spin polar- 
ization. If the interpretation of hydrogen-induced ferri- 
magnetism in HOPG [ll^ is correct, then with our tech- 
nique it should be possible to create the smallest con- 
trolled magnet ever possible by implanting hydrogen at 
the edges of carbon nanotubes. 

One of the authors (P.E.) thanks K. Kusakabe for cor- 
respondence. This work is supported by DFG Grant ES 
86/6-3. 
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